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We address the problems of fermions in light front QCD on a transverse lattice. We propose and 
numerically investigate different approaches of formulating fermions on the light front transverse 
lattice. In one approach we use forward and backward derivatives. There is no fermion doubling 
and the helicity flip term proportional to the fermion mass in the full light front QCD becomes an 
irrelevant term in the free field limit. In the second approach with symmetric derivative (which 
has been employed previously in the literature), doublers appear and their occurrence is due to 
the decoupling of even and odd lattice sites. We study their removal from the spectrum in two 
ways namely, light front staggered formulation and the Wilson fermion formulation. The numerical 
calculations in free field limit are carried out with both fixed and periodic boundary conditions on 
the transverse lattice and finite volume effects are studied. We find that an even-odd helicity flip 
symmetry on the light front transverse lattice is relevant for fermion doubling. 
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I. INTRODUCTION 

Light front Hamiltonian formulation of transverse lattice QCD |j], has many interesting features. With the gauge 
choice A + = A 0 + A 3 = 0 and the elimination of the constrained variable A~ = A 0 — A 3 , it uses minimal gauge 
degrees of freedom in a manifestly gauge invariant formulation exploiting the residual gauge symmetry in this gauge. 
So far encouraging results have been obtained in the pure gauge sector and in the meson sector with particle number 
truncation (for a recent review see, Ref. H). 

It is well known that fermions on the lattice pose challenging problems due to the doubling phenomenon. Light 
front formulation of field theory has its own peculiarities concerning fermions because of the presence of a constraint 
equation. As an example, the usual chiral transformation on the four component fermion field is incompatible with 
the constraint equation for nonzero fermion mass [|J . There have been previous studies of fermions on the transverse 
lattice Our approach in this work is quite extensive and aims to understand the origin of the doublers. We 

identify an even-odd helicity flip symmetry of the light front transverse lattice Hamiltonian, absence of which means 
removal of doublers in all the cases we have studied. This is closely related to the need to break chiral symmetry 
explicitly in the usual Euclidean formulation of lattice fermions. 

As we shall see later in this article, the presence of the constraint equation in light front field theory allows different 
methods to put fermions on a transverse lattice. It is worthwhile to study all the different methods in order to examine 
their strengths and weaknesses. In this work, we carry out a detailed numerical investigation of three methods in free 
field limit with special emphasis on finite volume effects. We also study the effects of imposing fixed and periodic 
boundary conditions which have significant effects in finite volumes. There are two important reasons to thoroughly 
study finite volume effects. First, for a reasonable size of Fock space, computing limitations will force us to be 
in a reasonably small volume when we deal with realistic problems. Second, the currently practiced version of the 
transverse lattice gauge theory uses linear link variables and recovering continuum physics is nontrivial. Finite volume 
studies are also important in this connection. 

In one of the approaches of treating fermions on the light front transverse lattice, we maintain as much transverse 
locality as possible on the lattice by using forward and backward derivatives without spoiling the hermiticity of the 
Hamiltonian. In this case doublers are not present and the helicity flip term proportional to the fermion mass in 
the full light front QCD becomes an irrelevant term in the free field limit. Thus in finite volume, depending on the 
boundary condition used, the two helicity states of the fermion may not be degenerate in the free field limit. However, 
we find that in the infinite volume limit the degeneracy is restored irrespective of the boundary condition. 
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In the second approach || , symmetric derivatives are used which results in a Hamiltonian with only next to nearest 
neighbor interaction when we take the free field limit. As a consequence even and odd lattice sites decouple and the 
fermions live independently of each other on the two sets of sites. As a result we get four species of fermions on a two 
dimensional lattice as excitations around zero transverse momentum (Note that this is quite different from what one 
gets in the conventional Euclidean lattice theory when one uses symmetric derivatives. In that case, doublers have 
at least one momentum component near the edge of the Brillouin zone.). The doublers can be removed in more than 
one way. We propose to use the staggered fermion formulation on the light front transverse lattice to eliminate two 
doublers and reinterpret the remaining two as two flavors. In this light front staggered fermion formulation, there 
is no flavor mixing in free field limit. But, in QCD, we get irrelevant flavor mixing terms. An alternative which 
removes doubling completely is to add the conventional Wilson term which generates many irrelevant interactions on 
the transverse lattice. Among them, the helicity flip interactions vanish but the helicity non flip interactions survive 
in the free field limit. 

The plan of this paper is as follows. Notation and conventions are presented in Sec. II. QCD Hamiltonian with 
forward-backward derivative is discussed and free field limit is studied in Sec. III. QCD Hamiltonian with symmetric 
derivative with its free field limit is considered in Sec. IV. Staggered formulation and reinterpretation of doublers are 
discussed in Sec. V. Removal of doublers via the Wilson term is studied in Sec. VI. We discuss the even-odd spin 
flip symmetry and its relation to the fermion doubling on the light front transverse lattice in Sec. VII. Finally Sec. 
VIII contains summary and conclusions. In appendix A we compare and contrast the forward-backward derivative in 
conventional lattice and light front transverse lattice for free fermion field theory. 


II. LIGHT FRONT PRELIMINARIES 


A. Notation and conventions 


The light front coordinates are x± = x° ±a; 3 , x ± = (x 1 , x 2 ), the partial derivative = 2^^, the gamma matrices 
^ = 7 ° ± y 3 and projection operators A ± = . x + is the light front time and x~ is the light front longitudinal 

coordinate. 

The Lagrangian density for the free fermion is 

£/ree = (2.1) 

Going to light front coordinates and using i/j^ = A^i/’, 

£/ree = V’ + ^'9 _ V’ + + — ip~\ia ± ■ 8 ± + 7 °m) , i/’ + — ip + \ia ± ■ 8 ± + 7 °m)ip~ . ( 2 . 2 ) 

One of the equations of motion from the above free Lagrangian is 

id + i\)~ = (ia ± ■ 8 ± + 7 °m)^ + , (2.3) 


which is actually a constraint equation because of absence of a time derivative. 
ip + is the dynamical fermion field and its equation of motion is given by 

id~ip + = ( ia 1 • d ± + 7°m)'i/’ _ - 


(2.4) 


One can remove from Eq. (2.4) using the constraint equation (Eq. (2.3)). The dynamical field ip + can essentially 
be represented by two components M such that 


ijj + (a 


-) = 


0 


-) 


(2.5) 


where 77 is a two component field. Its Fock expansion in the light front quantization with tranverse directions discretized 
on a two dimensional square lattice is given by 


V(x ,x) = 


Xx 


dk+ 


2(2ir)Vk + 


[b(k + , x, A)e 




+ d\k + , x, — A)e 2? 


( 2 . 6 ) 


where \x is ih e Pauli spinor, A = 1, 2 denotes two helicity states, x now denotes the transverse lattice points. 
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The canonical commutation relations are 

{b(k + , x, A), b\k' + , x', A')} = {d(k + ,x, X),cft(k ,+ ,x.', \')} 

= 2(2t r)k+ S(k+ - k ,+ ) 6\\> <5 XiX ' • (2.7) 


Using dk + ei k+ ( x v ) = 2(2tt)S(x —y ), we have, 

|? 7 (a; _ ,x), 77 t (|/ _ ,x')| = 1 <5 X>X / d(ar -y _ ) 


( 2 . 8 ) 


where 1 is a 2 by 2 unit matrix. 

We use Discretized Light Cone Quantization (DLCQ) [jIo| for the longitudinal dimension (—L < x~ < +L) and 
implement anti periodic boundary condition to avoid zero modes. Then, 


r)(x ,x) 


1 

V2L 


A 


53 Ml, x, \)e~ i '" lx ^ 2L ' > + rf(l, x, —X)e inlx 

(=1,3,5,... 


(2.9) 


with 


{fc(Z, x, A), 6 + (Z', x', A')} = {d(Z, x, A), dt(Z', x', A')} = S W 5 X ^ 6 X , X ,. (2.10) 

In DLCQ with antiperiodic boundary condition, it is usual to multiply the Hamiltonian P~ by so that H = j;P~ 
has the dimension of mass squared. 

In the following, for notational convenience we suppress x~ in the arguments of the fields. 


III. HAMILTONIAN WITH FORWARD AND BACKWARD DERIVATIVES 

A. Construction 

The fermionic part of the Lagrangian density is 

Cf = - m)ijj (3.1) 


with iD» =id fX - gA^. 

Moving to the light front coordinates, using the A + = 0 gauge and introducing the transverse lattice, 

Cf = i/j + \id~ — gA~)ip + + 

—a r D^ip + — a r D^.il>~ 

—TO7/) _ E°l/ ,+ ~ (3-2) 


Here r = 1,2 and D'l' b is the forward/backward covariant lattice derivative. Our goal here is to write the most 
local lattice derivative. That is why, instead of using the symmetric lattice derivative, in the above we have used the 
forward and backward lattice derivatives. However, the Hermiticity of the Lagrangian (Hamiltonian) requires that 
if one of the covariant lattice derivatives appearing in Eq. (3.2) is the forward derivative, the other has to be the 
backward derivative or vice versa. The covariant forward and backward derivatives on the lattice are defined as 


D( f ? 7 (x) = -[C/ r (x)? 7 (x + a r) - ry(x)] 


(3.3) 


and 

D b r r](x) = — [t]{x.) - f7 r f (x - ar)r](x - or)], (3.4) 

a 

where a is the lattice constant and r is unit vector in the direction r = 1,2 and Df = —D^. f7 r (x) is the group 
valued lattice gauge field with the property t/t(x) = t/_ r (x + or). Using the constraint equation 

id + ip~ = (ia r Df + 7 °m)ip + 1 


(3.5) 
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and finally going over to the two component fields g 


Cf = ip +1 (id — gA )ip + — tot/> + 7°^ — iip +1 a r D^.ip 
= — gA~)%l> + 

- V> +t [ia r D b r + 7°m] -7-— [ia s D{ + 7°m](/i + 

= (id~ — gA~)r) 

— [zcTr-Z?^ — ?TO]T^jr[*(T s -Df + im\rj . 


(3.6) 


(j 1 = < 7 2 and <72 = — <7i- So, we arrive at the Lagrangian density 


C f = g\id - gA )g - m 2 ^-^r/ 


( x )E^;m U r (x)g(x + ar) - g(x) 


mg' 


, id + . 


•- Y, [V( x + ar)U}. (x) — 77^ (> 


,11,, 

ar am v{x) 


- + -?7 t (x)]o- I .^<7 s [17 s (x)? ? (x + as) -?y(x)]. 


(3.7) 


In the free limit the fermionic part of the Hamiltonian becomes 


p fb ~ 


dx a 2 ^ H 


= dx 


T 


wV(x)—? 7 (x) 


- ^ f ( x ) Y r^( x + af ) _ 27 ?( x ) + ? ?( x - af ) 

r 

+ ^ + ( x ) Xl( arn ^)^:[i7( x + or) - 2 ? 7 (x) + t?(x- or)] 


(3.8) 


In order to get Eq. (3.8), we have assumed infinite transverse lattice and accordingly have used shifting of lattice 
points. The positive sign in front of the last term wou ld c hange if we had switched forward and backward derivatives. 

Because of the presence of <j r the last term of Eq. (3.8) couples fermions of opposite helicities. Note that it is also 
linear in mass. Such a helicity flip linear mass term is typical in continuum light front QCD. Here in free transverse 
lattice theory this term arises from the interference of the first order derivative term and the mass term, due to 
the constraint equation. This is in contrast to the conventional lattice (see Appendix A) where no helicity flip or 
chirality-mixing term arises in the free theory if we use forward and backward lattice derivatives. 

In DLCQ the Hamiltonian is given by, 


Hfb = H 0 + H h f 


(3.9) 


where 


H 0 = [ fot (^,z, A)6 (Z,z,A) +d f (Z,z, X)d(l,z, A)] 

z A l 

-EEEEEjxVxA 

z r A A' l 

b\l, z, A ')b(l, z + ar, A) — 2b\l, z, A ')b(l, z, A) + b\l : z, A ')b(l, z — of, A) 
+d" < (l , z, A r )d(l, z + of, A) — 2 d\l, z, A ')d(l, z, A) + d)(l, z, A ')d(l, z — ar, a) 


(3.10) 
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H hf = EEEEEy^'M^ 

z r A A' l 

b\l, z, A ')b(l, z + of, A) — 2&I(Z, z, X')b(l, z, A) + 6^(Z, z, X')b(l, z — or, A) 

+g^(Z, z, X')d(l, z 4- or, A) — 2d^(Z, z, X’)d(l, z, A) + d*(Z, z, X')d(l, z — or, tr) . 

(3.11) 


B. Absence of doubling 

Consider the Fourier transform in transverse space 


V(x ,x) = 


d 2 k 

W) 2 




where —— < k \, &2 < + -. Then the helicity nonflip part of Eq. (3.8) becomes 


(3.12) 


P nf = / dx~ 


Using 


d 2 k 


d 2 p 


)j4+M* )« 2 E ( 


(27r) 2 J (27r) 2 ^ kv ~ 'id+ 


i(k-p)-x 


m 


-yl 

„2 


gt l 




— 2 + e 


-*par 


: ^y(k-P)-x = (2 tt) 2 A 2 (k — p) 


(3.13) 


(3.14) 


we get, 


P nf = / <faT 


drk 


(27r) 2rk ^ '*0+' 


4(* )^T^k(a; ) 


E k r 


/ sin k r a/2 
\ k r a/2 


(3.15) 


where we have defined k r a = k • ra. Note that the sine function vanishes at the origin fci, &2 = 0 but does not vanish 
at the edges of the Brillouin zone k\. k -2 = ± -. 

Define k r = k r sll (. k a r ^/' 2 ■ In the naive continuum limit k r —> k r - 


Now, let us consider the full Hamiltonian (T8) including the helicity flip term. In the helicity space we have the 
following matrix structure for P + P~ (since P~ is inversely proportional to the total longitudinal momentum P + , we 
study the operator P + P~) 


™ 2 + ^Er sin 2 ¥ ^(*sin 2 ^ + sin 2 ^ 


— (isin 

a ' 

which leads to the eigenvalue equation 


2 k x a 


— sin' 


2 k y a 


m 


4E r sin : 

a z 


2 k r a 
2 


(3.16) 


. o 4 v-^ . 9 , 4m 

M 1 = to 2 4— 7 > sir —- ± —, 
a 2 2 a 


/E si 


k r a 


sin 


4 ■* , f 


(3.17) 


Third term in the above equation comes from the linear mass helicity flip term. If the mass m = 0, then it is 
obvious from Eq. (3.17) that A4 2 = 0 if and only if k\ = ^2 = 0. For nonzero m, one can also in general conclude 
that Ai 2 = m 2 only for the case k\ = k '2 =0. Thus there are no fermion doublers in this case (for physical masses 
am < 1). In the following for specific choices of momenta we elaborate on this further. 

If one component of the momentum vanishes, then 


4,1 


ka 


At 2 = m 2 H—(— ± m) sin 

a a 2 


(3.18) 
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where k is the non-vanishing momentum component. Thus for am = 1, irrespective of the value of k we get Af 2 = m 2 
which is unwanted. In general, for am > 1, Af 2 can become negative. It is important to recall that physical particles 
have m < ^ (the lattice cutoff) and hence are free from the species doubling on the lattice. With periodic boundary 
condition (discussed in the next subsection), allowed k values are k q = ±5n+T, 9 = 1,2, 3,..., n for 2n + 1 lattice 

sites in each direction. Let k\ = 0. For ma = 1.0, Eq. ( [3. If ) with the minus sign within the bracket gives M 2 = m 2 
for all values of k% and we get 2(2 n + l)-fold degenerate ground state with eigenvalue to 2 . 

The two spin states (spin up and down) are degenerate for k\ = /c 2 = 0. But if any one (or both) of the two transverse 
momenta is (are) nonzero then the degeneracy is broken on the lattice by the spin flip term proportional to to. So 
the total degeneracy of the lowest states for ma = 1.0 can be calculated in the following way: (a) k\ = /c 2 = 0 : 
Number of states =2(spin up and spin down), (b) k\ =0, ki^ 0 : Number of states =2 n and (c) k\ ^ 0, A '2 = 0 : 
Number of states =2n. Note that ki can have 2n nonzero values and there is no spin degeneracy for any nonzero ki. 
So, the total number of degenerate states = 2 + 2n + 2n = 2(2 n + 1). But if ma ^ 1 we cannot have to 2 eigen value 
for nonzero ki and we have only two (spin) degenerate states with eigenvalue m 2 . Again we see from Eq. ( 3.18 ) that 
if ma > 1, the kinetic energy term becomes negative and the eigenvalues go below to 2 . But ma > 1 means to > — 
(ultra violet lattice cutoff) and hence unphysical. 


C. Numerical Investigation 


We have investigated the effects of two types of boundary conditions: (1) fixed boundary condition and (2) periodic 
boundary condition. 


1. Fixed boundary condition 

For each transverse direction, we choose 2n + 1 lattice points ranging from — n to +n where fermions are allowed 
to hop. To implement fixed boundary condition we add two more points at the two ends and demand that the 
fermion remains fixed at these lattice points. Thus we consider 2 n + 3 lattice points. Let us denote the fermion 
wavefunction at the location s by u(s). We have u(s) ~ sin(s — 1 )ka with 7t(l) = u[2n + 3) = 0. Allowed values of k 
are (2 n + 2)k p a = pir with p = 1,2,3,...., 2n + 1 and k p = ( 2 w+ 2 )a -P- Thus the minimum k p allowed is ~ ( 2 ra+ 2 ) anc ^ 

maximum k p allowed is ~ [ 2 ^+ 2 ] • Tor example, for n = 1 we have ki = fc 2 = , etc. 


2. Periodic boundary condition 

Again, for each transverse direction, we choose 2n + 1 lattice points. We identify the (2 n + 2) th lattice point with 
the first lattice point. In this case we have the fermion wavefunction u(s) ~ e lska with the condition u(s) = u(s + L) 
where L = 2n + 1. Thus (2 n + 1 )k p a = ±27rp so that k p = ± ( 2 „+iv P> P = 0,1, 2,..., n. Thus the minimum k p allowed 
is 0 and the maximum k p allowed is 7 t 2 2 ” 1 . For n = 1, we have, ko = 0, fci = ±|7r, etc. 


3. Numerical results 


For the study of the fermion spectra on the transverse lattice, the longitudinal momentum plays a passive role and 
for the numerical studies we choose the dimensionless longitudinal momentum (l) to be unity which is kept fixed. For 
a given set of lattice points in the transverse space we diagonalize the Hamiltonian and compute both eigenvalues and 
eigenfunctions. _ 

First we discuss the results for H 0 given in Eq. ( |3.10 ). We diagonalize the Hamiltonian using basis states defined 
at each lattice point in a finite region in the transverse plane. Let us denote a general lattice point in the transverse 
plane by (a:,, ?/,:). For each choice of n (measure of the linear lattice size), we have —n < Xi,yi < +n. Thus for a 
given n, we have a (2n + 1) x (2 n + 1) dimensional matrix for the Hamiltonian. The boundary conditions do have 
significant effects at small volumes. For example, a zero transverse momentum fermion at finite n is/not allowed with 
periodic/fixed boundary condition. With fixed boundary condition, in the infinite volume limit, we expect the lowest 
eigenstate to be the zero transverse momentum fermion with the eigenvalue to 2 . In Fig. 1 we show the convergence 
of the lowest eigenvalue as a function of n towards the infinite volume limit in this case (to = 1 in Fig. 1). 
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For a zero transverse momentum fermion, the probability amplitude to be at any transverse location should be 
independent of the transverse location. Thus we expect the eigenfunction for such a particle to be a constant. At 
finite volume, with fixed boundary condition, we do get a nodeless wave function which nevertheless is not a constant 
since it carries some non-zero transverse momentum. All the excited states carry non-zero transverse momentum in 
the infinite volume limit. All of them have nodes characteristic of sine waves. The eigenfunctions corresponding to 
the first three eigenvalues are shown in Fig. 2 for the case of fixed boundary condition. With periodic boundary 
condition, for any n, we get a zero transverse momentum fermion with a flat wave function. 

Now, we consider the effect of hclicity flip term. With fixed boundary condition the lowest eigenstate has non 
vanishing transverse momentum in finite volume. In the absence of helicity flip term positive and negative helicity 
fermions are degenerate. The helicity flip term lifts the degeneracy. The splitting is larger for larger transverse 
momentum. In Fig. 3 we present the level splitting for the helicity up and down fermions as a function of n. As 
expected, the level splitting vanishes and we get exact degeneracy in the infinite volume limit. For the periodic 
boundary condition, the lowest state has exactly zero transverse momentum and we get two degenerate fermions for 
all n. 


IV. HAMILTONIAN WITH SYMMETRIC DERIVATIVE 
A. construction 


The symmetric derivative is defined by 


D r f/’ ± ( x ) = — [C/ r (x)'0 ± (x + ar) — t/_ I .(x)^/> ± (x — or)]. 


(4.1) 


In place of using forward and backward derivatives in Eq. (3.2), we use the above symmetric derivative for all lattice 
derivatives. Proceeding as in Sec. IIIA, we arrive at the fermionic part of the QCD Hamiltonian 


P sd = dx a 


^mV(x)—77(x) 

X 

- f dx~ a 2 ^2 < m^x) ^ ° r ^j+ [ c/ r(x)j?(x + ar) - U- r {x.)r](x - ar)] 

J x l, r 

— ?n—— ^2 [V( x — ar)<r r U r (x — ar) — r^(x + ar)a r U- r (x + ar)] tw—77(x) > 

Zi(Z tu 

r ) 

— J dx~ a 2 ^2 ~^~2 EE [ 7 I^( X — af)[/ r (x — ar) — (x + ar)t/_ r (x + ar) 


1 

i d + 


| U r (x)ri(x + ar) — [/_ r (x)?;(x — ar) 


In the free limit, the above Hamiltonian becomes 


P scL ~ 


dx a 2 


5Z| m V( x )^+^( x ) 

+ i “2 H irf ( x + «r) - (x - ar)] [r/(x + ar) - rj(x - ar )] j . 


(4.2) 


(4.3) 


In the free field limit the two linear mass terms cancel with each other. 

Using DLCQ for the longitudinal direction, we get 

Kd = ~H sd =- [H m + H k ] (4.4) 

7r 7r 

with 

H m = « 2 ™ 2 EEEj 

l a z 

[b'(l, z, cr)b(l 7 z, a) + d'(l, z, a)d(l, z, a)] 


(4.5) 





(4.6) 


H k 


EEEEy 

l a z r 


b'(l, z + ar, a)b{l, z + ar, <r) + &t(Z, z — ar, a)b(l, z — ar, a) 

— b'(l, z + ar, a)b(l, z — ar, a) — ZF(Z, z — ar, a)b(l, z + ar, a) 
+ d'(l, z + ar, cr)d(l, z + ar, a) + d'(l , z — ar, a)d(l , z — ar, a) 

— cP (/, z + ar, a)d(l , z — ar, a) — d^(Z, z — ar, a)d(l, z + ar, a) 


When we implement the constraint equation on the lattice and use symmetric definition of the lattice derivative, it 
is important to keep in mind that we have only next to nearest neighbor interactions. Thus a decoupling of even and 
odd lattice points occur and as a result we have two independent sub-lattices one connecting odd lattice points and 
the other connecting even lattice points. 

Let us now address the nature of the spectrum and the presence of doublers. 


B. Fermion doubling 


The Hamiltonian (fill) can be rewritten as 


P s d ~ 


dx a 2 Y [mV(x)^p??(x)-^ 2 a 2 [r? t (x)^ T E[ r ?( x + 2a? )]+ 7 ?( x “ 2a? )“ 2r ?( x )] 

x—even r 

+ Jdx~a 2 Y [TO 2 T? t ( x )^pr7( x )-^2a 2 [?7 t ( x )^j:E[ r ?( x + 2af )]+ 7 ?( x_2af ) _2? 7( x )] 

x—odd r 


(4.7) 


Clearly the Hamiltonian is divided into even and odd sub-lattices each with lattice constant 2 a. As a result, a 
momentum component in each sub-lattice is bounded by j- in magnitude. Again, going through the Fourier transform 
in each sub-lattice of the transverse space, we arrive at the free particle dispersion relation for the light front energy 
in each sector 





r 


k r a\. 


(4.8) 


For fixed k r , in the limit a —> 0 sin 2 k r a —> fc 2 and we get the continuum dispersion relation 


k 


k 


m 2 + k 2 
k+ 


(4.9) 


Because of the momentum bound of j- doublers cannot arise from ka = n. However, because of the decoupling of 
odd and even lattices, one can get two zero transverse momentum fermions one each from the two sub-lattices. Thus, 
for two transverse dimensions, we can get four zero transverse momentum fermions as follows: (1) even lattice points 
in x, even lattice points in y , (2) even lattice points in x, odd lattice points in y , (3) odd lattice points in x, even 
lattice points in y, and (4) odd lattice points in x, odd lattice points in y. Thus we expect a four fold degeneracy of 
zero transverse momentum fermions. 


C. Numerical Investigation 

1. Fixed boundary condition 

For each transverse direction, we have 2n + 1 lattice points where the fermions are allowed to hop. To implement 
the fixed boundary condition, we need to consider 2n + 5 lattice points. For one sub-lattice we have to fix particles 
at s = 1 and s = 2n + 5. We have, the wavefunction at location s, u s ~ sin (s — 1 )ka. We have, u s = 0 for s = 1. 
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We also need u s = 0 for s = 2n + 5. Thus (2 n + 4 )k p a = pn, with p = 1,2,3, n+ 1. For n = 1, allowed values of 

U orp U - _2I_ 27T 

e — 6a’ 6a * 

For the other sub-lattice, we fix the particles at s = 2 and s = 2n + 4. The wavefunction at location s, u s ~ 
sin (s — 2)fca. rt s = 0 for s = 2 and s = 2n + 4. Thus (2n + 2)k p a = pn with p = 1,2, 3, n. For n = 1, only allowed 
value of k is k = 


Combining the two sub-lattices, for n = 1, the allowed values of k are 


6a ’ 4a 


, and 


6a ’ 


2. Periodic boundary condition 

For a given n, fermions are allowed to hop at 2n + 1 lattice points in each transverse directions. Consider 2n + 3 
lattice points. For one sub-lattice (2 n + 3) rd lattice point is identified with the lattice point 1. For the other sub¬ 
lattice (2 n + 2) nd lattice point is identified with the lattice point 2. Wavefunction at point s, u s ~ e* s/co . We require 
e lka = e h 2n + 3 ) fea . Thus k p a = ± (2 n+ 2 ) > P = 0,1, 2,..., E±I. For n = 1, we have, ko = 0, ki = ±^. 

For the other sub-lattice we require e 2lka = e *( 2n + 2 ) feo _ Thus k p a = ±^p, p = 0,1, 2,..., For n = 1, allowed 
value of k = 0. Thus for n = 1, taking the two sub-lattices together, the allowed values of k are 0, 0, 


3. Numerical results 

The results of matrix diagonalization in the case of the symmetric derivative with fixed boundary condition are 
presented in Figs. 4, 5 and 6. In Fig. 4 we present the lowest four eigenvalues as a function of n. At finite volume, 
the four states do not appear exactly degenerate even though the even-odd and odd-even states are always degenerate 
because of the hypercubic (square) symmetry in the transverse plane. The four states become degenerate in the infinite 
volume limit. The eigenfunctions of the lowest four states are presented in Fig. 5 for n = 5. As they correspond to 
particle states, they are nodeless. All other states in the spectrum have one or more nodes. For example, in Fig. 6 
we show the eigenfunction corresponding to the fifth eigenvalue which clearly exhibits the node structure. 

With periodic boundary condition, for any n we get four degenerate eigenvalues corresponding to zero transverse 
momentum fermions. Corresponding wavefunctions are flat in transverse coordinate space. 


V. STAGGERED FERMION ON THE LIGHT FRONT TRANSVERSE LATTICE 


As we have seen in the previous section that the method of symmetric derivatives results in fermion doublers, we 
now consider two approaches to remove the doublers. In this section we study an approach similar to the staggered 
fermions in conventional lattice gauge theory. In the next section we will take up the case of Wilson fermions. 

In analogy with the Euclidean staggered formulation, define the spin diagonalization transformation 

V{x i,x 2 ) = {a 1 ) Xl {a 2 ) X 2 x(xi,x 2 ). (5.1) 

We see from the QCD Hamiltonian given in Eq. ( |4.2| ) with symmetric derivative that in the interacting theory 
(except for the linear mass term) and also in the free fermion limit, even and odd lattice sites are decoupled and 
the Hamiltonian is already spin diagonal. So, it is very natural to try staggered fermion formulation on the light 
front transverse lattice. In this section we shall follow the Kogut-Susskind formulation |n[ and present an elementary 
configuration spa ce analysis for two flavor interpretation. After the spin transformation the linear mass term in the 
Hamiltonian (4.2) becomes: 


f dx a 2 Y^ | m ^X f ( x )X^( x > r )^ [^(xMx + ar) - U- r (x)x(x- of)] 

x L r 

[X f ( x ~ ar)<fi(x,r)U r (x-ar) - x f (x + ar)0(x, r)U- r {x + a?)] j 


(5.2) 


where, </>(x,r) = 1 for r = 1 and <f>(x, r) = (— l) Xl for r = 2. After spin diagonalization, the full Hamiltonian in the 
free field limit becomes 


P sf = J dx |mV( x )^:x( x ) 
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+ 4^2 Y1 [x 1 ( x + ) - x f ( x - af)] [x( x + ar ) - x( x - af)] 

r 

1 -mx t (x)^ ^2 <f>(x, r)[x(x + af) - x( x - ar)] 


2 a 


+ ^ X f ( x - af)]^(x,r)^ r x( x )|- 


(5.3) 


The two linear mass terms cancel with each other in the free theory, but since they are present in the interacting 
theory we keep them to investigate the staggered fermions. 

Since all the terms in Eq. (pll) are spin diagonal, we c an put only a single component field at each transverse site. 
From now on, all the x’s and \ T ’s appearing in Eq. (5.3) can be taken as single component fermion fields. Thus we 
have thinned the fermionic degrees of freedom by half. Without loss of generality, we keep the helicity up component 
of x at each lattice point. 

Apart from the linear mass term in Eq. (^d|), all the other terms have the feature that fermion fields on the even 
and odd lattices do not mix. Let us denote (see Fig. 7) the even-even lattice points by 1, odd-odd lattice points by 
1', odd-even lattice points by 2 and even-odd lattice points by 2', and the corresponding fields by Xi etc. Then the 
first of the linear mass terms 


X f ( x ) yi+ ^ X ’ r M x + of ) ~ X( x “ af- )] 


id + 


(5.4) 


can be rewritten as (suppressing factors of a), 


Xi t ^p(ViX2 + V2X2O + X2 t ^(ViXi - V2X1') 
+X 1 ' t -^:(V i X2' - V2X2) + X2' f v^(ViXi' + V2X1) + B 


(5.5) 


where Vi and V 2 are the symmetric derivatives in the respective directions. Looking at Fig. 7 it is apparent that 
these Vi and V 2 can also be interpreted as a block derivative, i.e., finite differences between block variables. For 
example, V 1 X 1 = Xi(l>0) — Xi(0, 0). B represents the contribution from other blocks. 

Using Eq. (5.1), in terms of the nonvanishing components of 77 , we have 


m=Xi, r l2= 1X2, Vi' = *Xi'> W = -X 2 ' 


(5.6) 


An interesting feature of lattice points 1 and V is that fermion fiel ds m and 771 / have positive helicity. 772 and 772 ' have 
negative helicity. In terms of 77 fields the expression given in Eq. (|5.5|) can be written as 


r ?i t- 4 r(-*Vi 772 - V 2 772') + 77 ? 2 t T^-(Vi 77 i + iV 2 ? 7 i') 
io + io + 


1 


1 


+* r ?i' t_ Hz(~ v i ? 72' +*V 2 77 2 ) - t? 2 ' t T^-(-*V 1771 / +v 2 ?7i) +B. 


7 ( 9 + 


id + 


Now, 


77(1) — 77(0) = 77(?7(1) — 7?( 1)) + ^(77(1) + ??(—1) — 277(0)) 
= Vt 7(0) + iv 2 77(0) , 


(5.7) 


(5.8) 


77 ( 0 ) 7 ?( 1 ) = 77(77(1) 7 ?( 1 )) 77(77(1) T 7 ?( 1 ) 277 ( 0 )) 

= V77(0) - iv 2 77(0) 

where V and V 2 are respectively first order and second order block derivatives. So, we can write the expression 
as 


(5.9) 

d) 


-*(V 1? 72 - ^V 2 ?7 2 ) - (V 2 772' - ^772'} 


. t. 
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Let us introduce the fields 


+i?? 2 t T^j:{(Vi? 7 i + iVi 7 ?i) + *(V 2 7 ?i' - Jv^n)} 

+* 7 ?i' t ^+{ - (Vi??2' + iv?^)+i(V 2 »te + gV^)} 

-??2' t ^j:{ - i(Viryi/ - ^V???i<) + (V 2 ??i + ■ 


u 1 


u 2 


di 


d 2 


7f( r?2 + r?2 ') 

7! (r?2_r?2,) ' 


Then, the first order derivative term in Eq. ( 5.10 ) can be written as 

v) -rrrd r \7 r u + d) a r X7 r d = p —yr <r r V r / 


id + 

where, d = a 1 d and the flavor isospin doublet 


id + 


f = 


id + 


Similarly, we can write the second order block derivative term in expression ( 5.1(J| ) as 

ift — a 3 T r V 2 f 

2 J id+ rJ 

where, T r s are the matrices in the flavor space defined as 

1 = —ia , 1 = —ig . 

Similarly, the second term in Eq. ( |5.3D 

TW (x + ar) - p (x - or)] [ry(x + ar) - r](x - ar )] 

io + 


(5.10) 


(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 


+ v;/tTv;/ + j[v r /t + v’/t A. 


reads as 


The full Hamiltonian given in Eq. 

P~ f = f dx~a 2 


can now be written in two flavor notation as 




+ l i^rP-^a r T r Vlf + V?/t -L a r T r V r f)} 


i.c)} . 


(5.17) 


(5.18) 


The above simple exercise shows that applying the spin diagonalization on the symmetric derivative method, the 
number of doublers on the transverse lattice can be reduced from four to two which can be reinterpreted as two 
flavors. Although in the free case given by Eq. (5.18) the second and third lines are separately zero identically, we 
have kept these terms because in QCD similar terms will survive. These terms exhibit flavor mixing and also helicity 
flipping. The flavor mixing terms are always irrelevant. 
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VI. WILSON TERM ON THE LIGHT FRONT TRANSVERSE LATTICE 


Since doublers in the light front transverse lattice arise from the decoupling of even and odd lattice sites, a term 
that will couple these sites will remove the zero momentum doublers. However, conventional doublers now may arise 
from the edges of the Brillouin zone. A second derivative term couples the even and odd lattice sites and also removes 
the conventional doublers. Thus, the term originally proposed by Wilson to remove the doublers arising from ka = tt 
in the conventional lattice theory will do the job [[|. 

To remove doublers, add an irrelevant term to the Lagrangian density 


d£(x) = — V'( x )[^r(x)V’(x + ar) — 2i/’( x ) + t/_ r (x)r/;(x — a?)] 


where n is the Wilson parameter. This generates the following additional terms in the Hamiltonian ( 0 ): 


p- = 

x ^ r 

—4—— |V(x - ar)a r U r ht - ar) - rjUx + ar)dy[/_ I .(x + ar)] —r^x) > 

a 2 a L J io + I 

+ J dx~a 2 j ~T~ Y. Y [V( x ~ ar)I7 r (x — ar) + t?^(x + ar)?7- r (x + ar) 

x l. r s 

^U s (x)r](x + as) - U- S (x)r](x - as) 

— — — V |V^(x — ar)a r U r (x — ar) — 7^(x + ar)<5y£/_ r (x + ar) 

a 2 a L 

r s 

t/ s (x)? 7 (x + as) + U- S (pc)r](x - as) 

- f dx~a 2 ^ | A t ^ t ( x )^j: ^ V( X M X + ar) + [/_ r (x)??(x - ar)] 

^ x t r 

K x 1 ] 

+/i- ^2 V( x _ ar)U r {pi - ar) + ^(x + ar)J7_ r (x + ar)] tt^??(x) > 

/■ i K 2 r 

! 2J — ^ ^ ^(x - ar)t/ r (x - ar) + ^(x + ar)[/_ r (x + ar) 


av-^- [t/ r (x)r/(x + ar) - U- r (x)ri(x - ar)] 


r s 

l 

id + 


1 

id + 


|V s (x)? 7 (x + as) + U_ s (x)t 7 (x - as) 


In addition, the factor m 2 in the free term in (4.2) gets replaced by /i 2 = (m + 4^) 2 . 
In the free limit the resulting Hamiltonian goes over to 


( 6 . 1 ) 


( 6 . 2 ) 


P w = dx a 2 


^V( x )-^(x) 


+ a?) - 77 f (x - ar)]^ r ^-[77(x + ar) - r?(x - 


*9+ 2 a 


ar 


+ ^V (x + ar) - 2 ?y' 1 ' (x) + 77* (x - ar)] [77 (x + ar) - 2 ry(x) + ry(x - ar)] 

r 

-2—V (x)-i—[ry(x + ar) - 2?y(x) + ry(x - ar) 
a ' * 9 + 


(6.3) 


We rewrite the free Hamiltonian (6.3) as 


P w ~ Pd + Pqdi + Pqd2- 


( 6 . 4 ) 
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The diagonal terms are 

P d = J d x ~° 2 Y^ ^( x )^+ ? ?( x ) 

The nearest neighbor interaction is 

p om = - [ d ^ ffi2 EE 

J x r 

1 K 2 r 


\l 2 H-"b 8/i.Av— b 12 k; 2 — 

a z a cr 


(2/iK- + 4—) r 7 f (x)—ry(x + ar) + ?y t (x)—ry(x + ar) 


The next to nearest neighbor interaction is 


Pqd 2 — I dx a 


— J2 


E E { 4a 2 + a 2 } 


X r 
1 


^(x + ar)—— 77 (x — ar) + ryVx - ar)——rhx + of) 
?,£/+ '" n ' + 


id+ 


Using the Fourier transform in the transverse space, we get, 


P„„ = / dx 


d 2 k 








sin fc r 
k r a 


+2 ayLK^^k 


2 [ sinAva/2 
k r a /2 


+ a 2 K 2 E fe 


4 j sin k r a/2 
k r a/2 


(6.5) 


( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 


Note that, as anticipated, Wilson term removes the doublers because the lowest eigenvalue occurs only if all the k r 
are zero. 

In DLCQ, we have, 


Hd = [a 2 fi 2 + 1 + 8 a/iK + 12k; 2 ] ^ ^ ^ y 

l a z 

[&t(Z, Z, cr)6(Z, Z, cr) + d\l, Z, (T)rf(Z, Z, cr)], 


16.91 


and 


Horn = — [2fta/r + 4K 2 ] EEEEy 

l cr z r 

b'(l, z, cr)6(Z, z + ar, a) + b'(l, z, a)b(l, z — ar, cr) 

+ d'(l, z, cr)rf(Z, z + ar, cr) + cF(Z, z, cr)cZ(Z, z — ar, cr) 


( 6 . 10 ) 


p od 2 - -« 2 ]EEEEy 

l cr z r 

b'(l, z + ar, a)b(l, z — ar, cr) + Zr (Z, z ~ ar, a)b(l, z + ar, cr) 

+ gF (Z, z + ar, (t)cZ(Z, z — ar, a) + d'(l, z — ar, a)d(l, z + ar, cr) 


(6.11) 
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A. Numerical Investigation 

1. Boundary condition 


With the Wilson term added, we do not have decoupled sub-lattices. We have both nearest neighbor and next- 
to-nearest neighbor interactions. Since with fixed boundary condition, the lowest four eigenvalues are not exactly 
degenerate in finite volume, it is difficult to investigate the removal of degeneracy by the addition of Wilson term. 
With periodic boundary condition, for a lattice with 2n + 1 lattice points in each transverse direction, we identify 
the (2 n + 2) th lattice site with the first lattice site. Then for the Hamiltonian matrix we get the following additional 
contributions. 


( 



NN 0 
V N NN 


NN N \ 
0 NN 



( 6 . 12 ) 


The matrix elements NN = — j + k 2 and N = —2 apn — 4k 2 . For a given n, the allowed values of k are k p = ± , 

p = 0,1, 2,. Thus for n = 3, we expect multiples of =§• apart from 0. For n = 5, apart from 0, allowed values of k 

are multiples of 4y. 


2. Numerical results 

Since the Wilson term connects even and odd lattices, the extra fermions that appear at zero transverse momentum 
are removed once Wilson term is added as we now have nearest and next to nearest neighbor interactions. For 
large n, we get the expected spectra but, numerical results suggest that the finite volume effect is larger for small k 
which is obvious because k is a mass-like parameter. For example, with periodic boundary condition, for n = 3, for 
k = 1.0,0.5,0.4, we get the expected harmonics but not for k = 0.1. The situation is similar for n = 5. For n = 10, 
expected harmonics emerge even for n = 0.1 but not for k = 0 . 01 . 


VII. DOUBLING AND SYMMETRIES ON THE LIGHT FRONT TRANSVERSE LATTICE 

Because of the constraint equation which is inconsistent with the equal time chiral transformation in the presence 
of massive fermions, we should distinguish between chiral symmetry in the equal time formalism and in the light front 
formalism. For example, the free massive light front Lagrangian involving only the dynamical degrees of freedom is 
invariant under 75 transformation. On the light front, helicity takes over the notion of chirality even in presence of 
fermion mass which can be understood in the following way. 

In the two component representation J 9 ] in the light front formalism, let us look at the objects ^ and i/j^. We 
have 


t> + W=(’ ,( o ) ) (7.D 

with 

The projection operators are Pr= 1 + y 5 ) and Pl = \{ 1 — 7 ' 5 ) with 



(7.3) 
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Then 


and 


1>t = P R ^ + = 



= Pl^ + = 



(7.4) 


(7.5) 


Thus i/Jr = PriP + represents a positive helicity fermion and = PliP + represents a negative helicity fermion, even 
when the fermion is massive. This makes sense since chirality is helicity even for a massive fermion in front form. 
This is again to be contrasted with the instant form. In that case the right handed and left handed fields defined by 
ip R = P R if) = |(1 +7 5 )ip and ipL = Pl^P = ^(1 — J 5 )ip contain both positive helicity and negative helicity states. 
Only in the massless limit or in the infinite momentum limit, ip R becomes the positive helicity state and ipL becomes 
the negative helicity state. 

As a passing remark, we would like to mention that in continuum light front QCD there is a linear mass term that 
allows for helicity flip interaction. 

In lattice gauge theory in the Euclidean or equal time formalism, because of reasons connected to anomalies (the 
standard ABJ anomaly in vector-like gauge theories), there has to be explicit chiral symmetry breaking in the kinetic 
part of the action or Hamiltonian. Translated to the light front transverse lattice formalism, this would then require 
helicity flip in the kinetic part. A careful observation of all the above methods that get rid of fermion doublers on the 
light front transverse lattice reveals that this is indeed true. 

In particular, we draw attention to the even-odd helicity flip transformation 


v{xi,x 2 ) -> {a 1 ) Xl (a 2 )°° 2 r](x 1 ,X 2 ) 


(7.6) 


that was used in Sec. [v] for spin diagonalization. It should also be clear that the form of the above transformation is 
not unique in the sense that one could exchange dy and 02 and their expo nent s x\ and X 2 could be changed by ± 1 . 

Note that the Hamiltonians Pj b given in Eq. (3.8) and P~ given in Eq. ( |6.3| ) that do not exhibit fermion doubling 
are not invariant under the transformation Eq. (7.6). On the other hand the Hamiltonian P~ d given by Eq. ( |4.3| ) that 
exhibits fermion doubling is invariant under this transformation. 


VIII. SUMMARY AND CONCLUSIONS 

The presence of the constraint equation for fermions on the light front gives rise to interesting possibilities of 
formulating fermions on a transverse lattice. We have studied in detail the transverse lattice Hamiltonians resulting 
from different approaches. 

In the first approach, forward and backward derivatives are used respectively for ip + and ip~ (or vice versa) so that 
the resulting Hamiltonian is Hermitian. There is no fermion doubling. The helicity flip (chiral symmetry breaking) 
term proportional to the fermion mass in the full light front QCD becomes an irrelevant term in the free field limit. 
With periodic boundary condition one can get the helicity up and helicity down fermions to be degenerate for any 
transverse lattice size n. With fixed boundary condition, there is a splitting between the two states at any n but the 
splitting vanishes in the large volume limit. 

In the second approach, symmetric derivatives are used for both ip + and %p~. This results in four fermion species. 
This is a consequence of the fact that the resulting free Hamiltonian has only next to nearest neighbor interactions 
and as a result even and odd lattice sites get decoupled. One way to remove doublers is to reinterpret them as flavors 
using staggered fermion formulation on the light front. In QCD Hamiltonian, it generates irrelevant flavor mixing 
interactions. However, in the free field limit, there is no flavor mixing. Another way to remove the doublers is to add 
a Wilson term which generates many extra terms in the Hamiltonian. In the free field limit, only the helicity nonflip 
terms survive. The Wilson term couples even and odd sites and removes the doublers. Numerically, we found that in 
small lattice volumes it is preferable to have not too small values of the Wilson mass n/a. 

We have tried to understand the fermion doubling in terms of the symmetries of the transverse lattice Hamiltonians. 
We are aware that there are rigorous theorems and anomaly arguments in the conventional lattice gauge theories 
regarding presence of fermion doublers. In standard lattice gauge theory, some chiral symmetry needs to be broken 
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in the kinetic part of the action to avoid the doublers. On the light front, chirality means helicity. For example, a 
standard Wilson term which is not invariant under chiral transformations in the conventional lattice gauge theory, 
is chirally invariant on the light front in the free field limit. The question is then why the Wilson term removes the 
doublers on the light front transverse lattice. The argument that there is nonlocality in the longitudinal direction 
cannot hold because, in the first place, having nonlocality is not a guarantee for removing doublers and secondly, 
there is no nonlocality on the transverse lattice. One, therefore needs to find a reasoning that involves the helicity 
in some way. We have identified an even-odd helicity flip symmetry of the light front transverse lattice Hamiltonian, 
absence of which means removal of doublers in all the cases we have studied. 

Our interest also lies in studying finite volume effects on a transverse lattice. As we have emphasized, there are 
important issues to be understood since (a) any realistic Fock space truncation will force us to work with relatively 
small volumes because of the limited availability of computing resources and (b) the currently available transverse 
lattice formulation uses linear link variables and recovering continuum limit is nontrivial. We have investigated the 
effects of fixed and periodic boundary conditions, which are significant in finite volumes. 

Among the many possible extensions of this work, it will be interesting to study the various QCD Hamiltonians 
and to compare the resulting spectra. 
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APPENDIX A: FORWARD-BACKWARD DERIVATIVE IN CONVENTIONAL LATTICE THEORY 


In this appendix we follow Ref. ||l2| . In discretizing the Dirac action in conventional lattice theory the use of forward 
or backward derivative for d M leads to non-hermitian action. The hermiticity can be preserved in the following way. 
In the chiral representation 


1 

o 


' 0 a 1 ' 

5 

r i o 

-I 0 

, Y = 

-a 1 0 

, 7 = 

1 

o 


The Dirac operator in Minkowski space 




0 -ia^dy 

0 


(Al) 


(A2) 


where, = (I,<r), a 11 = (/, — a). For massive Dirac fermions, this leads to the structure 


- icr^dnipR - mipL 

(A3) 

—ia^dnipL - mtl’R- 

(A4) 

For discretization we replace d M in Eq. (|A3|) by forward derivative 

Aj£ = (8y tX +fj, — Sy, x )/a 

(A5) 

and in Eq. (A4) by backward derivative 

A^ = (Sy tX — 5y tX -n)/a. 

(A6) 

This leads to the structure 

m = i7 M A* - *7 M 75A“ - m 

(A7) 

which results in hermitian action. Here, 

A m = Yy,x+y — Sy tX -ij.)/2a 

A m = (Sy^x+y T fiy,x—y. 2dy ;3 ,)/2tt. 

(AS) 
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Note that irrelevant helicity nonflip second order derivative term is produced in this method of discretization. In 
contrast, the corresponding term in the transverse lattice depends linearly on m and flips helicity. One can trace this 
difference to the presence of the constraint equation in the light front theory. 
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FIG. 2: Eigenfunctions of first three states for the case of no doubling. n=5 
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FIG. 5: Eigenfunctions of first four (degenerate) states for the case of fermion doubling. n=5 
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FIG. 6: Eigenfunction corresponding to the fifth state. n=5 















